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We study self-adjoint operators defined by factorizing second order differential 
operators in first order ones. We discuss examples where such factorizations intro- 
duce singular interactions into simple quantum mechanical models like the harmonic 
oscillator or the free particle on the circle. The generalization of these examples to 
the many-body case yields quantum models of distinguishable and interacting par- 
\ tides in one dimensions which can be solved explicitly and by simple means. Our 

considerations lead us to a simple method to construct exactly solvable quantum 
many-body systems of Calogero-Sutherland type. 
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1. Introduction. 



In this paper we define and study peculiar self-adjoint extensions of simple differential 
operators like the one denning the quantum mechanical harmonic oscillator. We also find 
the explicit solution of some of the resulting quantum mechanical models and discuss 
their physical interpretations. These examples provide amusing illustrations that the 
mathematical theory of Schrodinger operators jKj RSl| IRS2} IBS . and in particular the 



distinction between formally hermitian and self-adjoint operators, can have important 
physical consequences even in elementary quantum mechanics. Moreover, these examples 
allow a straightforward generalization to the case of an arbitrary number of particles, 
and the resulting many-body models describe strongly interacting systems which can be 
solved by elementary means. These models include, as special cases, the ones studied 
by Girardeau [G], which have received considerable interest recently since it has become 
possible to realize and study them experimentally [F] . We stress that these models can be 
solved even in the general case of distinguishable particles, and that they can be studied 
by much simpler methods than other well-known integrable quantum many-body systems 
|LL| IY| IC2| |S] . Moreover, our considerations lead us to a remarkably simple and general 
approach to exactly solvable quantum many-body systems of Calogero-Sutherland type 
|C2t El IQPj which seems new. 
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The starting point for this work was a paradox for the differential operator —d 2 + x 2 
defining quantum harmonic oscillator, — oo < x < oo and d x := d/dx. The paradox arises 
if one confronts the formal factorization 



and Q*Q > 0, with the well-known fact that the groundstate eigenvalue of the harmonic 
oscillator Hamiltonian H = —d 2 + x 2 is 1. It is interesting to note that a variant of this 
paradox was already discussed by Calogero in 1969; see [01] . second paragraph after Eq. 



The paradox above is only one of many, and they are resolved by observing that 
each such factorization gives rise to a different quantum mechanical Hamiltonian defined 
by the Friedrichs extension. These examples have in common that a simple differential 
operator with an obvious factorization also allows for alternative factorizations where 
singularities are introduced, and this implies a restriction on the domain of functions on 
which the resulting self-adjoint operator is defined. This restriction has a natural physical 
interpretation as strong repulsive interaction. 

We will present several examples leading to interesting quantum mechanical model 
with additional interactions and which are exactly solvable in the sense that all eigen- 
functions and eigenvalues can be found explicitly and by simple means. 1 In particular, for 
each non-negative integer n we will find a quantum mechanical Hamiltonian H n which is 
defined by the differential operator of the harmonic oscillator differential operator — d 2 +x 2 
but is bounded from below by 2n + 1. Moreover, in addition to the standard harmonic 
oscillator H = H also the Hamiltonian Hi defined by the factorization in Eq. ([T} is ex- 
actly solvable, while the other H n > 2 are only partially solvable, i.e., only a finite number 
of eigenstates and eigenvalues can be found exactly. Two other examples are provided 
by the differential operator —d 2 — 7(7 + l)cosh _2 (x) often referred to as Poschl-Teller 
potential, and the differential operator — d 2 for — ir < x < ir describing a free particle on 
the circle. In the latter example we will show that, for every a > 0, the factorization 



defines a unique self-adjoint operator H a bounded from below by a 2 . Moreover, this 
Hamiltonian is exactly solvable if and only if a is integer or half integer. 

While we believe that these examples are of pedagogical interest, our first main result 
is that much of what we say can be naturally extended to the many-body case. This 
yields many-body systems describing distinguishable particles moving in one dimension 
and interacting through a strong local two-body interactions, and also these models can 
be solved exactly and by simple means. 

Our considerations described above lead us, somewhat unexpectedly, to our second 
main result which is on Calogero- Sutherland type systems: As we will see, all our one- 

1 To illustrate what we mean by this we recall that the spectrum of the Hamiltonian defined by the 
differential operator —d% on all C°°-functions on an interval and general boundary conditions is given 
by a transcendental equation in general which can be solved only numerically. These examples are not 
exactly solvable in this sense. 





(3.24). 



d 2 = (d x + a cot(ax)) (— d x + a cot(ax)) + a 2 



2 



particle examples mentioned above are related to Schrodinger operators H = —d x + V(x) 
which have exact eigenf unctions of the form 

Mx) = Pn{<n{x))z~ w{x \ n = 0,1,2,..., N max , (2) 

where P n is a polynomial of order n, rj and w are real- valued C 2 functions, and N max < oo. 
These quantum mechanical models are exactly solvable in the sense that these eigenfunc- 
tions ip n and the corresponding eigenvalues E n can be found explicitly. While we will 
restrict our discussions to our three examples: harmonic oscillator, Poschl- Teller poten- 
tial, and free particle on the circle, it is important to note that there are many more: in 
Tabled we list the well-known examples associated with the classical orthogonal polyno- 
mials, but even this list is not complete. Our second main result is a surprisingly simple 
argument showing that, to any such model, the corresponding quantum-many body given 
by the Schrodinger operator 

N 

H N = J2{- 9 l+ v ^)+ x ( x -^ E w &> x k) ( 3 ) 

3=1 l<j<k<N 

and the following particular two-body interaction 

_ j{xf + j{yf 2j(x)w'(x) - r]"{x) - 2j{y)w'{y) + rfjjj) _ 

[,y) (v(x)-v(y)) 2 riW-viv) [) 

defines an exactly solvable quantum many-body system; r]'(x) := d x r](x), and Wo is an 
arbitrary constant which can be fixed to some convenient value. More specifically, this 
arguments leads to a simple proof of the following: 



Proposition: Let H = —d x + V(x) be a self-adjoint Schrodinger operator on L 2 ([a,b}), 
— oo<a<b< oo, with exact eigenf unctions as specified in Eq. ff. Then the cor- 
responding many-body Schrodinger operator Hn in Eqs. (E^-j^P defines a quantum-many 
body systems with the following exact ground state, 

V (x 1 ,...,x N )=e-^ w ^ H (vix^-vix^y, (5) 

l<j<k<N 

provided N < N max + 1 and A > is such that \l/o is square-integrable. The corresponding 
eigenvalue is determined by the one-particle eigenvalues E n as follows, 

N-l 

£ = A 2 J2( E n - Eo) + NE - -A(A - l)N(N - 1)W . (6) 

n=0 



(A concise proof can be found in the last paragraph of Section 6.) 



In fact, we found that it is possible compute all other bound states and corresponding 
eigenvalues for all these systems exactly and in a unified manner |HLj . As we will discuss 
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Name 


1pn(%) 


^(xj E n 


T 
1 


Jacobi 


^SlIl^J [COS 2 J * n [COSX) 


4a 2_l 4/3 2 ~l / . a+/3+l\^ 
16sin*§ 1 16 cos 2 f ^ 2 ' 


TT 




1 Olll Jo 1 ^ XL \ wO J> 1 


a{a ~ l) in 1 nl 2 


III 


Chebyshev 


T n (cosx) 


n 2 


IV 


Laguerre 


e- 2 /2 x ( a +i)/2 L W( x2) 


4n + 2a + 2 x 2 + 



Table 1: Solutions of the Schrodinger equation (— <9 2 + V(x) — E n )ip n (x) = associated 
with classical orthogonal polynomials; adapted from Table 22.6 in |ASj . 



in Section 7, special cases of this result yield various well-known exactly solvable systems 
of Calogero-Sutherland type |C2t El IQPt IBFj . but there are other interesting examples 
which we have not seen discussed in the physics literature. 

The plan of this paper is as follows. In the next section we resolve in detail the paradox 
described above using the theory of Friedrichs extensions. This provides a general recipe 
for finding other such examples. In Section 3 we present our one-body examples, and in 
Section 4 we discuss their physical interpretation. In Section 5 we discuss the extension 
of our results to the many-body case. Our method to find quantum-many body systems 
with exact groundstates is explained in Section 6. We end with a few remarks in Section 7. 

2. Factorizations and selfadjoint extensions. 

We first recall some standard facts from the operator theory in Hilbert spaces. Let A be 
a symmetric operator in a separable Hilbert space TC with scalar product (•, •) (linear in 
the second argument) and norm || ■ ||, and let T>(Aq) be the domain of the operator Aq 
which is a dense subset in 7i. We assume that Aq is bounded from below, i.e., there exists 
a constant E EM. such that 

(A u,u) > E \\u\\ 2 , ueV(Ao). 

The sesquilinear form ao[u, v] = (A u,v) defined for u,v G d[a ] = V(A ) can be com- 
pleted with respect to the norm defined by 

a [u,u] + {E Q + l)\\u || 2 . (7) 

Such a completion is denoted by a and its domain by d[a\. By standard arguments (see 
for example fBSJ) there is a unique selfadjoint operator A corresponding to the closed 
form a whose domain T>(A) C d[a\. This operator is called the Friedrichs extension of the 
operator Aq. 

Let V be a real function defined on open subset Q C M. N , and let if be a self-adjoint 
Schrodinger operator in L 2 (Q) with some boundary conditions. 

Hu(x) = -Au(x) + V(x)u(x). (8) 

This operators generates a quadratic form h[u, u] defined on functions from d[h]. We 
assume that H has discrete spectrum {E n }^ =Q , with the corresponding orthonormal in 
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L 2 (Q) system of eigenfunctions {u n }^ =0 where K < oo (i.e. H may also have continuous 
spectrum). It can be easily checked directly that the differential expression —A + V(x) 
can formally be factorized as 

- A + V = Q* n Qn + E n , (9) 

where 

Qn = -V, + (10) 

It is well-known that the ground state «o can be chosen strictly positive and the factoriza- 
tion @ does not give any problem. However, considering higher eigenfunctions we obtain 
singularities in the right hand side of (fTUj) due to the zero set S of u n , S = {x : u n (x) = 0}. 
In order to define the selfadjoint operator corresponding to the operator Q*Q we begin 
with the quadratic form 



h n \u,u] 



-V x u + - — -u 



Vu„ 2 N 
ax, 



whose domain d[h n ] consists of functions u G d[h] such that u(x) = for x £ S. 

The corresponding Friedrichs extension gives us a selfadjoint operator which we denote 
by H n . Although 

H n = Q* n Q n + E n = -A + V (11) 

coincides with H = — A+V defined by (JHJ) regarding its action on functions as a differential 
expression, its domain is very different and consists of functions not only just satisfying 
the boundary conditions of the original operator H but being equal to zero on S. This 
automatically makes the domain d[h n ] smaller compared with the domain of the original 
operator H and "lifts up" the spectrum. 

In particular, this answers the question concerning the formal inequality which follows 
from (fTTj) . that is 

H = Q* n Q n + E n > E n . 

The equality in the right hand side of this latter equation is wrong because the domains 
of the operator H and Q* n Q n are different and therefore H ^ Q* n Q n + E n as Hilbert space 
operators. 

3. Some examples. 

Here we shall give some simple examples where the above construction can be worked out 
more explicitly. 

3. a. Harmonic oscillator. 

Let H be the operator of the harmonic oscillator on M, 

H = -d 2 x +x 2 . (12) 

The eigenvalues and eigenfunctions of this operator are well-known and equal to E n = 
2n + 1, n = 0, 1, . . . , and 

M$ = H n {x)e- x2 l\ (13) 
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where Ti n are Hermite polynomials, 

H n (x) = (2 n n\)- 1/ \- 1/i (-l) n e x2 d2e- x \ 
As in fill)) the operators H n can be defined via 

H n = Q* n Q n + E n = [d x + |^) (-d x + |^) + En- 

According to Section 2 the latter expression coincides with —d 2 + differential 
expression but not as a selfadjoint operator. In particular, when n = 0, 

— = -x, 

and we obtain the standard factorization of the harmonic oscillator 

h = -dl + x 2 = (d x - x) (-a, - x) + i. 

If n = 1 then 

— = -ac + -, 

ipi x 

and we obtain the factorization 

H x = -dl + x 2 = (d x - x + (-^ - x + ~) + 3. 

It is interesting to note that not only Hq but also Hi defines an exactly solvable 
system in the sense that all its eigenfunctions can be computed explicitly: Using that the 
eigenfunctions ip n (x) of H in Eq. (|13|) ff. are even and odd for even and odd integers n, 
respectively, we conclude that the eigenfunctions u n (x) of Hi are given by 

u 2n +i(x) = ip 2 n+i(x), u 2n {x) = 4>2n+i(x) sign(x), u = 0, 1, 2 . . . (14) 

with the usual sign function sign(x): the odd eigenfunctions of Hi are the same as for 
the standard harmonic oscillator, but the even ones are very different. The spectrum of 
Hi thus has multiplicity two, E 2n = E 2n +i = 4n + 3. 

The other exotic harmonic oscillator Hamiltonians H n > 2 are only partially solvable: as 
discussed in Section 4 below, we can determine the (n + l)-fold degenerate ground state 
of H n > 2 explicitly, but it is not possible to compute its excited states analytically. 

3.b. Poschl-Teller potential. 
Let 



H = ~dl - (15) 



7(7 + 1) 
cosh 2 x 

for 7 > 0. It is well-known that this operator has [7] (= largest integer smaller or 
equal to 7) square integrable eigenfunctions ip ni n = 0, 1, . . . , [7] — 1, with corresponding 
eigenvalues 

E n = -{ 1 -nf. (16) 
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These eigenfunctions are given by 

v r_iv 

4 s (n - 2s)! (7 + s - n)\s\ 



ipJx) = cosh" 7 (a;) V (-1Y-, ^1 — rnsmh n - 2s (x) (17) 



0<s<n/2 

with r! := T(r + 1) for non-integers r, as can be checked by straightforward computations. 

In particular, the lowest eigenvalue of this operator is E = — 7 2 and the corresponding 
eigenfunction is ipo(x) = cosh~ 7 (:r). Then 

ip' 

— = — 7tanh(x) 
Wo 

and therefore one can obtain the following factorization 

#0 = -d 2 x - 2i2±i2 = f dx _ 7 tanh(a:) N ) (-d x - jtznh(x)) - 7 2 
cosh x V /V / 

defining the standard selfadjoint operator Hq = H given by the differential operator 
— d 2 +7(7 + 1) cosh _2 (x). The second eigenvalue is E\ = —(7 — l) 2 and the corresponding 
eigenfunction is given by ipi(x) = cosh~ 7 (:r) sinh(x), for 7 > 1. Then 



= coth(x) — 7tanh(x) 



which implies that 

H x = -d 2 x - = (q x + coth(x) - 7 tanh(x)) 

cosh x \ / 



x (-d x + coth(x) - 7tanh(x)j - (7 - i f 



which defines another selfadjoint operator Hi, and similarly for H n >2- Again it is possible 
to compute all eigenfunctions and eigenvalues of Hi in terms of the ones of H , similarly 
as in the harmonic oscillator case. 



3.c. Free particle on the circle. 

This example is slightly more subtle. For arbitrary a > we define a self-adjoint operator 
on L 2 ([— 7r, 7r]) by the factorization 

H a = — 9 2 = (d x + acot(ax))(—d x + acot(aa;)) + a 2 

which is obtained with the odd eigenfunctions sin(aa;) of — d 2 . The subtle point now 
concerns the correct boundary conditions at x — ±7r. According to our general theory 
in Section 2, we need to start with the self-adjoint operator H = —d 2 with boundary 
conditions such that sin(o;a;) is indeed its eigenfunction, and it is not difficult to see that 
these conditions are as follows, 

u'{— 7r) = u'(tt) and u(n) — u(— 7r) = — 2 tan(mr)it'(7r) (18) 
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for functions u in the domain of H. Thus only for integers a this corresponds to periodic 
boundary conditions, and only in this case can we interpret H as Hamiltonian of a particle 
on a circle. We conclude that the eigenfunctions of H a are defined by the same boundary 
conditions at x — ±7r together with the requirements that they vanish at all points x 
where sin (ax) = 0. 

If we interpret H Q as Hamiltonian of a particle moving on a circle, then we not only 
have infinitely strong repulsive interactions at the points x where sin(ax) vanishes but, 
in addition, a singular interaction at the point x = ±7r whose strength depends on a so 
that it vanishes at integer values of a and becomes infinitely strong if a is half integer. 

In particular, if a = 1/2, then the operator 

#1/2 = ~dl ={d x + \ cot(x/2)) (-d x + | cot(z/2)) + \ 

is defined on functions u satisfying a Dirichlet boundary condition at zero (from the 
singularity of cot at zero) and Neumann boundary conditions at — 7r and tt, 

u (0) = u'(-tt) = u'(tt) = 0. 

Our discussion above implies that, for integers and half-integers a, one can determine 
all eigenfunctions and eigenvalues of H a explicitly and by elementary computations, but 
otherwise part of the spectrum is determined by a transcendental equation: the eigen- 
functions are determined by the conditions u(kn/a) = for integer k and the boundary 
conditions in Eq. (fTSj) . and the latter are non-trivial unless 2a is an integer. 

For example, for a = 1 these eigenfunctions are given by 

u± tTl (x) = 6(±x) sin(na;) 

with n = 1,2,.... The corresponding eigenvalues are n 2 and two fold-degenerate. For 
a = 1/2 the explicit eigenfunctions are 

u±,„ = Q(±x) sin((n + l/2)ar/2) n = 0, 1, 2, . . . 

with corresponding two-fold degenerate eigenvalues (n + 1/2) 2 . 

We finally mention that it is amusing to visualize the a-dependence of the singular 
interactions and the spectrum of the Hamiltonian H a . 

4. Physical Interpretation. 

We first give a physics argument providing an intuitive resolution of the paradox men- 
tioned after Eq. (JTJ). 

The paradox arises due to the singularity in the factorization in Eq. (|T|) at x = 0. One 
can regularize this singularity and consider the following family of Hamiltonians, 
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with a regularization parameter e > 0. By a straightforward computation one finds 

e 2 (l + 2e 2 + 2x 2 ) 



x 2 + e 



2\2 



Even though the second term vanishes for e — 0, it is obviously relevant in the limit 
£ | 0: Recalling that S E (x) = e/(x 2 + e 2 )7i is a sequence of C°° functions converging 
to the delta distribution as e j 0, we see that this term gives an additional singu- 
lar interaction of infinite strength at x = 0: formally, the second term converges to 
lim £ j (^^(x)) 2 = ii 2 5(0)5(x) = lim g ^ +00 g 6(x). This interaction forces the eigenfunc- 
tions of the Hamiltonian lim £ jo to vanish at x — 0. It is plausible that this latter 
Hamiltonian is identical with the one defined by the factorization in (fT|). and this is indeed 
the discussed below. 

As shown in the previous section, the eigenfunctions of Hamiltonians H = —d 2 + V(x) 
on L 2 ([a, b}), —oo<a<b<oo, defined by the factorization 

H n = -d 2 x + V(x) = (d x + u' n (x)/u n (x)) (-d x + u' n {x)/u n {x)) + E n 

can be fully characterized by the (simple) zeros x n j of the eigenfunction u n of — d x + V(x): 
u n (x n ,j) = for j = 1, 2, . . . , n implies that u(x n j) = for all eigenfunctions u of H n . 
One thus can interpret H n as the strong coupling limit g — > oo limit of the following 
Hamiltonian, 

n 

= -d 2 x + V(x) + Y,9S(x-x n ,). 

Indeed, the eigenfunctions u of this latter Hamiltonian are defined by —u" + Vu = Eu 
for x x n j, together with the boundary conditions 

u(x nJ + 0) = u(x n>j - 0), u'(x nJ + 0) - u'(x nJ - 0) = gu(x nij + 0), 

and for g — > oo the latter reduce to u(x n j) = 0. Thus the factorization corresponds 
to adding infinitely strong repelling delta-potentials at the locations of the zeros of the 
eigenfunction u n . This implies that the ground state of H n is n + 1-fold degenerate: all 
states uoj(x) which are equal to u n (x) in an interval between two adjacent zeros x n j and 
zero otherwise are eigenfunctions of H n with eigenvalue E n . However, the excited states 
cannot be found in such a simple manner in general. 

5. Exactly solvable many-body systems in ID. 

We now show that the examples above allow a straightforward extension to the iV-particle 
case, and this gives rise to exactly solvable interacting many body systems of distinguish- 
able particles. 

5. a. Interacting particles in a harmonic oscillator potential. 
Our first example is for the operator 

N 

H N = J2(- d l+ x *) ( 19 ) 

3=1 
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on R . This self-adjoint operator is defined by the standard factorization 

N N 

h n = + xj) = £ (d Xj - Xj ) {-d Xj - x 3 ) + N, 

3=1 3=1 

and it describes an arbitrary number, N, of noninteracting, identical particles moving in 
an external harmonic oscillator potential. The eigenfunctions of obviously are given 
by products of the one particle eigenfunctions defined in Section 3. a, 

N 

*Pn(x) = Y[^ n] (x j ), Uj = 0,1,2,... , 

3=1 

where we use the shorthand notation x = {xx, . . . , x^) and similarly for n, and the 
eigenvalue is 

N 

£„ = 5^(2^ + 1). (20) 

3=1 

Using the huge degeneracy of these eigenfunctions we can construct eigenfunctions of 
which vanish at all coinciding points Xj = Xk by antisymmetrization, 

N 

f n (x) = (-^Tl^njixpj), < m < n 2 < . . . < n N < oo. (21) 

Pes N j=i 

These are exactly the fermion eigenfunctions of Hn, of course: exchanging any two par- 
ticles gives a minus sign. The fermion wave function with the lowest eigenvalue is 

fo(x) = const det (%_i(^)) (22) 

l<j,k<N 

where const is some irrelevant non-zero constant which we will later fix to some convenient 
value. To compute /q we use the explicit form of the harmonic oscillator eigenfunctions 
in Eq. 0, 

f no (x)=e~^ x V 2 det (7^-1(2*)). 

Since the Hermite polynomials are of the form Ti n (x) = c n x n + . . . with dots some linear 
combinations of Ti, m<n and c n non-zero constants, we can replace in this determinant 
TCj-i(xk) by cj^ixjr , and a well-known identity of the Vandermonde determinant implies 
that 

fo{x) = e~ E ^ X V 2 Yl ( Xj - Xk 



3 

l<j<k<N 

AT— 1 , 
3= 



is the fermion ground state of the operator H N with the eigenvalue ^ ?=0 (2j + 1) = N 2 

"j 

N 



(we fixed const so as to absorb factors of Cj and a possible sign). Then obviously 

1 d 



-f {x) = + 



f (x)dxj ^[ x j- x k 
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We thus conclude that the differential operator defining Hn also can be factorized as 
follows, 

N N N 

H NA = Y,(-d 2 X] +x 2 J ) = Y,Q* J Q 3 +N 2 , Q, -<),, - .r, ■ £ . (23) 

3=1 j=l k=l X i Xk 

k^j 

and this defines another self-adjoint operator. 

We now discuss how to determine all eigenf unctions f(x) of Hn,i] see jO] for a similar 
discussion. As in Section 2 we conclude that the eigenfunction have to vanish at coinciding 
points, 

f(xi, X2, ■ ■ ■ , %n) = if Xj = Xk for some j ^ k, j, k = 1, 2, . . . , N, (24) 

and thus all fermion eigenfunctions /„, of Hn in Eq. (J2T|) are also eigenfunctions of Hn,i 
with the same eigenvalues. It is not difficult to see that all these are C°° eigenfunctions of 
.ffjy,i) but they obviously do not span the full Hilbert space L 2 (WL N ) but only its fermion 
subspace. To find the missing eigenfunctions we divide the space R in fundamental 
domains which can be parameterized by permutations P G Sn, 

A P := {a? = (xi, . . . , Xn) £ K | — oo < Xp\ < Xp 2 < ■ ■ ■ < Xp^ < oo} . (25) 

There are N\ distinct such domains, and the closure of the union of all these domains 
equals to full space, Upgs^Ap = W N . Moreover, the boundaries of these domains corre- 
spond to the hyperplanes Xj = Xk- We thus need to find all functions which are C°° and 
obey the eigenvalue equation for in the interior of all these domains Ap and vanish at 
all their boundaries. The functions in Eq. (}2*T]) obey both these conditions, but since the 
derivatives of eigenfunctions need not be continuous at coinciding point also the functions 

rPf \ f fn(x) for X G Ap n „ no 

ax) = otherwig / , < n, < n 2 < . . . < n N < oo, P G S N (26) 

are legitimate eigenfunctions: they obey the correct eigenvalue equation with the eigen- 
value in Eq. (J2*U|) at non-coinciding points, and they also obey all boundary conditions in 
Eq. ()24j). It is not difficult to verify that these eigenfunctions f n {x) span the full Hilbert 
space L 2 (R N ). 

Similarly as discussed in Section 4, we can interpret H^ i as the strong coupling limit, 
g — > oo, of the following Hamiltonian 

N 

3=1 l<j<k<N 

describing particles moving on the real line in the harmonic oscillator potentials and 
interacting with two-body delta function potentials. This corresponds to the fact that 
the quadratic form of the operator Hn,i is split into iV! independent problems in L 2 (Ap) 
with Dirichlet boundary conditions. 
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5.b. Interacting particles in a Poschl-Teller potential. 



We now consider the operator 

N 

^ = £(-^-7(7+1)^) (27) 



describing N noninteracting particles in a weakly confining cosh 2 -potential. For 7 > N 



the fermion ground state of this Hamiltonian is fo as in (j22|) but now with the eigenfunc- 
tions ip n in (|17|) . and by an argument as in Section 5. a we conclude that 

N 

fo(x) = JJcoslf ] I (sinh(x i ) - sinh(x fc )) 

j=l l<j<k<N 



2 



with corresponding eigenvalue 

£ o = ~ $^(7 - nf = ~j(2N -1)(N-1) + jN(N - 1) - N>y 

n=0 

We thus can define the following Hamiltonian different from Hn, 

N 

H N>1 = J2Q*3Qi + £o (28) 
i=i 

with 

Qj = ~d Xj + U X J = -d Xj - 7 tanh(x J ) + ■ ?° S ^ ] (29) 
Jo fcli sm H x j) ~ sm(xfe) 

As in Section 5. a we can interpret Hn,i as the strong coupling limit of the Hamiltonian 
Hn in ()27|) with two body delta interactions enforcing that the eigenf unctions vanish 
at coinciding points Xj = Xk, and all eigenvalues and eigenfunctions of this interacting 
many-body Hamiltonian can be determined explicitly from the fermion wave functions 
and corresponding eigenvalues of the non- interacting Hamiltonian H^. The example here 
is, however, quite interesting from a physical point of view since, in addition to many- 
body bound states, there are also scattering states and, in general, one has mixed states 
where only some of the particles are trapped and the others are free. 

5.c. Interacting particles on the circle. 
Let 

N 

H N = -A = -J2d 2 X] (30) 
i=i 

on dimensional torus T^. The right hand side of (|3T)|) can be considered as a trivial 
factorization of the periodic Laplace operator. We note that 



^„(x) = ( e lx ) n e- i{N - 1)x/2 , n = 0, ±1, ±2, 
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all are eigenfunctions of the corresponding one-particle Hamiltonian — <9 2 with eigenvalue 
(n — ^(N — l)) 2 , and thus fo as defined in (J22j) is an fermion eigenstate of the operator 
(JHUJ) on T N with the eigenvalue - |(JV- 1)) 2 = ±N(N 2 - 1). By straightforward 

computations we find that 

fo{x)= \ \ sml(xj-Xk), 

l<j<k<N 

and therefore, by analogy with (|23|). we can consider the factorization 

N N 

0*0. + 

12 



N N 

Hs A = - E = E + t^ n ( n2 - 1) ( 31 ) 



i=i i=l 



where 



1 ^ 
<2i = - -rd Xj f = -d Xj + V I cot - z fc ). (32) 

/o *=i 

Mi 

Note that the singularities enforce that eigenfunctions of H^ i vanish at points Xj = x^ 
modulo 2tt, and thus we can interpret H^ i as describing particles moving on a circle of 
length 27r, — 7r < Xj < 7T, interacting with a strong two-body delta interaction. 

As in the previous examples this factorization splits the operator H^,i into an orthog- 
onal sum of operators in fundamental domains parameterized by permutations P G Sn, 
and the eigenfunctions of -Hat,i can all be determined explicitly as antisymmetrized prod- 
ucts of the eigenfunctions of the corresponding non-interacting Hamiltonian H^. 

6. A method to find Calogero-Sutherland type systems. 

As already mentioned, our considerations in the previous section are closely related to 
arguments which are well-known in the context of Calogero-Sutherland type systems |Clj . 
In this sections we explore this relation in more detail. Through that we are naturally lead 
to a surprisingly simple approach to a large class of Calogero-Sutherland type systems. 

One well-known method to find quantum integrable systems is to write down a many-body 
state of Jastrow form 

N 

*ow=n e_, " w n 9(xj,x k )\ a>o, 

j=l l<j<k<N 

for some C 2 functions w and g, and then construct a many-particle Hamiltonian 

- 1 

h n = E Q i q j - Q * = ~ d *s + *r d *i ^° 

with S some real, conveniently chosen constant. Then obviously defines self-adjoint 
operator > £ , and Qj^o = Vj implies that ^ is an exact groundstate of H N with 
eigenvalue £ . In general, such a Hamiltonian H N describes a system of particles with two- 



13 



and three-body interactions, and only for very particular functions w and g the three-body 
terms cancel. These cases are of particular interest since they lead to physically interesting 
Hamiltonians of the form as in Eq. (J3J) which, in many cases, are exactly solvable. Famous 
such examples are w(x) = \x 2 ) g(x,y) = (x — y) and w(x) = 0, g(x,y) = sin|(x — y) 
corresponding to the models discovered by Calogero |C2| : 



^ = ^(-^+* 2 )+A(A-l 
and Sutherland S\: 



j=l l<j<k<N v y > 



N 



H N = -Y,d 2 Xj +\{\-l) 



1 



^ 2sin 2 i(x 7 - 

'i<k<N 2 ^ J 



i=x i<j<fc<iv" 01 " 2V " J Xj) 

respectively. Olshanetsky and Perelomov found interesting variants of the Calogero- 
Sutherland models related to root systems; see |OPj and references therein. 



Obviously, the Calogero- and Sutherland models are very closely related to our ex- 
amples discussed in Sections 5. a and 5.c, respectively. It is interesting to note that our 
considerations there provide an alternative proof of the key identity leading to these mod- 
els. Indeed, if we check the identity in (J2~3"j) by a direct computation, we find that it is 
equivalent to the well-known functional identity which cancels the three-body terms in 
the Calogero model: 

_ + cycl = 

[Xj - X k ){Xj - Xi) 

for j, k,£ all different, where L +cycV is short for adding the two terms of the same form 
obtained by cyclic permutations of the indices, (j, k,£) — > (k,£,j) and (j, k, £) — > (£,j, k). 
However, since we have a direct proof of (|2*Hj) based on a basic identity of the Vandermonde 
determinant, we found an alternative proof of this functional identity. Similarly, the 
identity in Eqs. (|31j) - (|32j) is equivalent to the well-known functional identity 

cot(\(xj — Xk)) cot \{xj — xi) + cycl = —1, 

used to cancel the three-body terms in the Sutherland model. 

It therefore is natural to ask: is there a Calogero-Sutherland type system corresponding 
to our example in Section 5.b? We now can use a similar method as above to find such a 
system: we can use Eqs. (J2BlH(I2El) to derive the functional identity 

cosh 2 (a; 1 ) 

+ cycl = l, 



(sinh(x-,) — sinh(a;fc)) (sinh(x-,) — sinh(a^)) 
which then can be used to cancel three-body terms and prove that 

N 

^o(x) = Y\ cosh~ 7 (:rj) Y\ (smh(xj) — smh(xk)) X (33) 

j=l l<j<k<N 
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is the exact ground state of the Hamiltonian 



H N = y (-81 - ^±in + \(X - 1) y oosh^x^ + cosh 2 ^) 
3 cosh {xj)J 1 J^< N (sinh( Xj ) -sinh(x fe )) 

Moreover, it is straightforward to show that all other bound states and corresponding 
eigenvalues can be found using Sutherland's method [Sj- 2 

It is interesting to note that our argument can be easily extended to a rather large 
class of examples: assume a one-particle Hamiltonian H = —d 2 + V(x) which has exact 
eigenstates of the form 

V>„(z) = P n (v(x))e- W{x \ n = 0,l,2,... 

with eigenvalue E n , where P n is a polynomial of order n and w and r] are C 2 functions. 
Then f in (J22j) is obviously an eigenstate of the many-body Hamiltonian 

»< + V( Xj )) 
j 

with eigenvalue £ = J2f=o E n . As in Section 5 we can conclude that 

N 

j=l l<j<k<N 

and using this formula for / the identity 

N 1 

+ v i?i)) = E ^ + ^ = + 7^/0 

i=i i=i /o 

implies a functional identity expressing three-body terms through two-body terms. We 
now compute 

h = J2q*q j , Qj = -d Xj + ±d Xj % 

with \l/o in ©, and using this functional identity we can cancel all three-body terms and 
thus find that Tt = H N — £ with H N in Q-® and £ in (JEj). Since obviously H^o = 
and H. > 0, we obtain the result stated in the Proposition in Section 1. Table Ogives 
other interesting examples of this result. We plan to give a comprehensive discussion of 
all these models and their solution in a future publication HL . 

7. Final comments. 

We should mention that the model in (J34j) is closely related to the well-known CV-variant 
of the Sutherland model defined by 

J sum { X j)J i<f^<N^ 2sm 2^-^') 

1 



2 E. Langmann, unpublished. 



2 sin =(xj + Xj) 
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see jOPj and references therein. This Hamiltonian describes identical particles in an 
external, repulsive potential and with repulsive two-body interactions, and it only has 
scattering states (purely continuous spectrum). One can show that by shifting Xj — > 
Xj + 'vk/2, this Hamiltonian is converted to the one in (|34j) . up to an additive constant. 
However, these shifts changes the physics of the model dramatically: the Hamiltonian 
in describes interacting particles in a weakly confining potential, and it thus has, 
in addition to bound states and scattering states, also mixed states where some of the 
particles are confined and others free. To our knowledge, this situation has not been 
studied in the context of Calogero- Sutherland type systems. We therefore think it would 
be interesting to study this model in detail which, however, is beyond the scope of the 
present paper. 

There are other special cases which seem less known in the physics literature even though, 
to our opinion, they also of interest in physics. The examples listed in Table [T] lead to 
systems which were previously known in the mathematics literature; see |BF[ IDXj and 
references therein. We should mention that these latter systems are closely related to 
the BCn variants of the Calogero- (IV) and Sutherland model (I III) |OPj . respectively. 
However, our result nevertheless seems useful since it is the starting point of an alternative, 
systematic and unified approach to these systems and their solutions, including examples 
which are hardly known [HLJ. 
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